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ABSTRACT 

We test the luminosity function of Milky Way satellites as a constrain for the nature of 
Dark Matter particles. We perform dissipationless high-resolution iV-body simulations 
of the evolution of Galaxy-sized halo in the standard Cold Dark Matter (CDM) model 
and in four Warm Dark Matter (WDM) scenarios, with a different choice for the 
WDM particle mass (m w ). We then combine the results of the numerical simulations 
with semi-analytic models for galaxy formation, to infer the properties of the satellite 
population. Quite surprisingly we find that even WDM models with relatively low 
m w values (2-5 keV) are able to reproduce the observed abundance of ultra faint 
(M v < —9) dwarf galaxies, as well as the observed relation between Luminosity and 
mass within 300 pc. Our results suggest a lower limit of 1 keV for thermal warm dark 
matter, in broad agreement with previous results from other astrophysical observations 
like Lyman-a forest and gravitational lensing. 
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1 INTRODUCTION 

The inflationary cold dark matter scenario gives a clear pre- 
diction for the initial fluctuation spectrum responsible for 
the formation of the Large Scale Structure in the Universe, 
with considerable power down to very small scales. As a con- 
sequence we expect the mass function of dark matter haloes 
to rise steeply towards low masses. N-body simulations have 
indeed revealed that in Cold Dark Matter (CDM) models, 
all Galaxy-sized haloes {Mdm ~ 1O 12 M0) should contain a 
large number of embedded subhaloes that survive the col- 
lapse and virialization of the parent structure (e.g Springel 
et al. 2008). Although the predicted number of substruc- 
tures is in reasonable agreement with observed luminosity 
functions (LFs) in cluster sized haloes, in Milky Way (MW) 
sized haloes the number of predicted sub-haloes exceeds the 
number of observed satellites by at least an order of mag- 
nitude (Klypin et al. 1999, Moore et al. 1999). This tension 
between models and observations represents, together with 
the detailed reproduction of the star formation histories (e.g. 
Salvadori et al. 2008) and mass profiles of local dwarf galax- 
ies (e.g. Walker et al. 2009), one of the most relevant ques- 
tions for the present day theories of galaxy formation and 
evolution. 
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A possible cosmological solution to this discrepancy is 
to replace cold dark matter with a warm species (WDM, 
Bode, Ostriker & Turok 2001 and references therein). The 
warm component acts to reduce the small-scale power, re- 
sulting in fewer galactic subhaloes and lower halo central 
densities (Colin et al. 2000,2008). One possible WDM candi- 
date is a sterile neutrino which exhibit a significant primor- 
dial velocity distribution and thus damp primordial inho- 
mogeneities on small scales (e.g. Hansen et al. 2002, Abaza- 
jian & Koushiappas 2006 Boyarsky, Ruchayskiy, Shaposh- 
nikov 2009). Since ground based experiments may be far 
from directly studying in detail the nature and properties 
of the actual DM particle, it is of fundamental importance 
to determine astrophysical constraints on the maximum free 
streaming length of any warm candidate. 

Limits on the mass of dark matter particles can be ob- 
tained from several astrophysical observations: theoretical 
phase space density studies combined with observations of 
stellar dynamics in Milly Way satellites (e.g. Boyanovsky, de 
Vega, Sanchez 2008; de Vega Sanchez & Sanchez 2009), lu- 
minosity function of high redshift QSOs (Song & Lee 2009), 
abundance of dwarf galaxies in the Local Volume (Zavala 
et al. 2009) and the size of (min)-voids around the Local 
Group (Tikhovov et al. 2009). Perhaps one of the most pow- 
erful tool for constraining the matter power spectrum are 
Lyman-a forest observations (neutral hydrogen absorption 
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in the spectra of distant quasars, Narayanan et al. 2000, Viel 
et al. 2005). Lyman-a observations allows the possibility to 
studying the power spectrum down to small scales and over 
a large range of redshifts. Using HIRES data, a lower limit of 
m w — 1.2 keV has been reported by Viel et al. (2008). Even 
tighter limits can be obtained combining Lyman-a observa- 
tions with SDSS results (see Boyarsky et al. 2009a and ref- 
erences therein) . These results have been recently confirmed 
by Miranda & Maccio (2007) in an independent estimation 
of lower limits for the mass of the WDM particle based on 
QSO lensing observations. 

In the last five years our knowledge on the number 
and properties of satellite galaxies around our the MW has 
tremendously increased thanks to the results from the Sloan 
Digital Sky Survey (SDSS). The homogeneous sky coverage 
of the SDSS provided the first determination of the volume 
corrected MW satellite luminosity function down to lumi- 
nosities as faint as 100 Lq (Koposov et al. 2008). Exploiting 
this advancement in the observational knowledge of our own 
Galaxy, several theoretical works have revised the problem 
of satellite number density, coming to the conclusion that it 
is possible to reconcile the observational evidences with the 
predictions of the standard (L)CDM scenario (e.g. Koposov 
et al. 2009, Maccio et al. 2009, Li et al. 2009a). 

The aim of the present work is to expand the results 
presented in Maccio et al. 2010 (M10, hereafter), trying to 
use a combination of observational data and theoretical pre- 
dictions to infer significant constraints on the allowed mass 
of any warm dark matter particle. We start from the findings 
of M10 and we re-simulate one of the DM haloes presented 
in that work, in a WDM scenario for different choices of m w . 
We then combine the results with the best fit Semi Analyt- 
ical Model (SAM) codes, as defined in M10, to compare the 
resulting prediction for the luminosity function of MW satel- 
lite in order to infer a lower limit on the WDM candidate 
mass. 



2 SIMULATIONS AND SAMS 

In this paper we combine merger trees extracted from very 
high resolution N-body simulations, describing the hierarchi- 
cal assembly of a MW-like halo, with SAM techniques, to 
predict the relationship between the dark matter (sub)haloes 
and observable galaxy properties, allowing us to make a di- 
rect and detailed comparison with observational data from 
SDSS. We perform LCDM simulations and analyze them (in- 
cluding the definition of DM haloes and the reconstruction 
of their detailed merger tree) using the same tools described 
in full detail in M10 and we refer the reader to that paper 
for a more detailed discussion. In particular, for this work 
we select one specific halo (namely the Gl halo in M10) and 
we resimulate it for a suite of WDM models with particle 
masses m w — 10, 5, 2, 1 keV. The dark matter particle mass 
is nid — 4.16 x 10 5 /i _1 M©, and gravitational softening of 355 
h~ x pc, with w 3 x 10 6 particles within the virial radius. 

To generate initial conditions for WDM, we define a 
rescaled version of the CDM power spectrum using a fit- 
ting function that approximates the transfer function asso- 
ciated to the free streaming effect of WDM particles (Viel 
et al. 2005). We do not include the effect of a non zero pri- 
mordial velocity dispersion for WDM particles in our simu- 



lations, this because even in our more extreme WDM model 
(1 keV) the rms of the random velocity component at the 
starting redshift of the simulations (z « 40) is of the order of 
2 km sec -1 . This is much smaller than the typical velocities 
induced by the gravitational potential itself. Hence we do 
not expect that neglecting this random velocity component 
could alter our results. In the following, we will not show 
the results for the m w = 10 keV model since it turned out 
to be indistinguishable from LCDM. 

In order to predict the expected luminosities of satel- 
lite galaxies, we combine the results of the iV-body simu- 
lations with two state-of-the-art SAMs, namely the Kang 
et al. (2005, 2006, see also Kang 2008) model and MOR- 
GANA(Monaco et al. 2007, updated in Lo Faro et al. 2009). 
In M10 we studied in full detail the effect of the various 
physical processes in shaping the luminosity function of MW 
satellites. We determined that suppression of gas infall by 
a photo-ionizing background, supernova feedback and tidal 
destruction are the most relevant processes responsible for 
the agreement between theoretical predictions and data. We 
characterize our best fit models as follows, (i) We regulate 
star formation efficiency in low mass haloes by shutting off 
gas cooling in structures with virial temperature below 10 4 K 
(due to the inefficiency of H2 cooling), (ii) We suppress 
hot gas accretion in low mass haloes according to photo- 
ionization background, reionization and filtering mass argu- 
ments (Kravtsov et al. 2004). (iii) Stellar feedback is mod- 
eled as a function of halo circular velocity. For the purposes 
of this work we keep the parameters fixed at the same values 
in the best fit models of MloQ The two SAMs provide con- 
sistent prediction with respect to the dependence of the LF 
on the WDM particle mass, and lead to similar conclusions. 
For a sake of simplicity, in the following of the paper we will 
focus and show only results referring to the Kang model. 



3 RESULTS 

The effects of WDM are clearly visible in figure [T] where we 
show the dark matter density map within 360 kpc for the 
2 extreme models: LCDM (upper panel) and LWDM with 
m w — 1 keV (lower panel). The difference in the number of 
substructures can be quantified by looking at the cumulative 
satellite mass function which is shown in figure [2] In order 
to identify bound subhaloes we use the Amjfl halo finder 
(Knollmann & Knebe 2009). 

In the WDM scenario not only the total number but 
also the formation history of dark matter structures is ex- 
pected to be different with respect to LCDM. Due to the 
lack of power on small scales (sub) haloes will both form and 
be accreted onto the main halo at different times. This ef- 
fect could leave a footprint on the luminosity of the satellites 
hosted by the subhaloes. Subhaloes will have a different mass 
at the time of reionization and this will affect the extent of 
gas suppression. Moreover, their accretion redshift onto the 
main halo will be different, with relevant implications for the 

1 A more detailed and critic discussion of how the variation in 
SAM parameters may affect the final luminosity function can be 
found in M10. 

2 The Amiga Halo finder (ahf) can be freely downloaded from 
http: / /www. popia.ft.uam.es/AMIGA 
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Figure 1. Dark matter density map within a sphere of radius H = 360 
kpc. LCDM and LWDM results (m w — 1 keV) are shown in the upper and 
lower panel respectively. 



star formation history of the satellites (in the SAMs we con- 
sider there is no new gas accretion onto subhaloes). Figure 
[3] shows the distribution of accretion (z acc ) and formation 
(zform) redshifts, where this latter quantity is defined at 
the time when the halo virial temperature exceeds 10 4 K 
(allowing gas cooling), for the LCDM and LWDM (m m — 2 
keV) models. As expected, gas cooling and star formation 
start later in the WDM scenario, with almost no haloes with 
Zform = 11 (upper panel). At the same time, they seem to 
have (on average) a later accretion time in WDM models 
with respect to LCDM (lower panel). Those two effects act 
on opposite ways in shaping the stellar content of the galaxy 
hosted by the dark matter (sub) halo, and, as we will see 
next, the net effect on the satellite luminosity function is 
almost negligible. 

Our main results on the MW satellite luminosity func- 
tion are presented in figure [4] In each panel, we compare 
the original results for the Gl halo in the LCDM cosmology 
(shown as a solid line) with the predictions corresponding to 
the different LWDM cosmologies (dotted lines). We decide 



to show only Gl in this work since, among the four haloes 
presented in M10, it was the one giving the best agreement 
with the observed satellite luminosity function. This makes 
the comparison between the observed and simulated LF in 
WDM models more straightforward. It is evident from the 
figure that the number density of satellites of different lu- 
minosity is reproduced for a WDM particle m w > 2 keV as 
well as for the standard cosmology. Only when considering 
m w — 1 keV we detect a significant discrepancy with the ob- 
servational data, with only nine predicted satellites fainter 
than My — —9. This result can be used to set a lower limit 
on on the warm dark matter particle mass (m m > 1 keV), in 
broad agreement with the results coming from the analysis 
of the galaxy power spectrum. 

It is worth noting the peculiar behavior of the bright 
end of the luminosity function. There is no satellite brighter 
than My — — 12 in the m w = 5 keV and m w = 2 keV real- 
izations, while the predictions for the m w = 1 keV model are 
in broad agreement with the observational data. We checked 
that this peculiar behavior is not related to the modeling of 
the dynamical friction or tidal stripping in the SAMs, but 
it is related to the different properties of Gl merger tree 
in the various LWDM cosmologies. As an effect of the re- 
duced small-scale power, the statistical properties of DM 
haloes change: haloes are on average less concentrated (Eke 
et al. 2001) and both the accretion time and the mass of a 
dark matter halo at that time are modified with respect to 
LCDM. Given the small number statistics associated to the 
bright end of the satellite LF for a galaxy-sized halo (less 
than 4 objects in total), this result is not totally unexpected. 
For this reason we stress that the stronger constraints come 
from the comparison of the faint end, where the number 
statistics are high enough to reveal any significant modifica- 
tion of the LF with respect to the LCDM result. 
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Figure 3. Lower panel: Distribution of accretion redshifts z acc for sub- 
haloes with M{z acc ) > 10 9 M Q . Solid (black) 1 inc shows results for the 
LCDM model, while the dot-dashed (red) line is for the LWDM with 
m w — 2 keV. Upper panel same as lower panel but for the formation red- 



shift z 



defined as the redshift at which the halo virial temperature 



exceeds 10 K. 



Mateo et al. (1993) and Strigari et al. (2008, see also 
Walker et al. 2009) have noted that all satellites with lu- 
minosity between 10 3 and 1Q 7 Lq have a common mass of 
~ 1O 7 M within a radius of 300pc (Mo. 3). Moreover, within 
this inner region, all objects result to be dark matter domi- 
nated. Maccio, Kang & Moore (2009, see also Li et al. 2009b, 
Okamoto & Frenk 2009) showed that this relation is not to- 
tally unexpected in a LCDM Universe. Given the expected 
modification in the properties of DM haloes as a function of 
redshift in a WDM model, it is interesting to check if our 
simulated satellites are still able to reproduce the observed 
normalization and small scatter of the M0.3 — L relation. 

Satellite luminosities are direct outputs of the SAM, 
while for computing Mo. 3 we used the same approach de- 
tailed in Maccio, Kang & Moore (2009). At time of accre- 
tion of each satellite we compute the parameters (r s and 
S c ) that describe its density profile assumed to be NFW 
(Navarro et al. 1997, see Maccio et al. 2008 for more de- 
tails). The present value of Mo. 3 is then computed by inte- 
grating the density profile, under the assumption that r s and 
5 C do not evolve with time. Figure [5] shows the distribution 
of the the mass within 300 pc as a function of the satellite 
luminosity. The upper panel presents results for a WDM 
with m m = 5 keV (blue open squares), which turn out to 
be indistinguishable from the LCDM ones (red points) and 
hence in good agreement with the observational findings of 
Strigari et al. 2008. The lower panel shows results for our 
most extreme WDM model (m w = 1 keV): in this case sim- 
ulated satellites suggest a stronger correlation between their 
inner mass and luminosity which is in (slight) disagreement 
with the flat distribution of the observational data. Never- 



Figure 4. Milky Way satellites luminosity function in LCDM and LWDM 
models. In each panel solid line refers to the prediction for LCDM cosmol- 
ogy, while dotted lines refer to the predictions for LWDM realizations, as 
indicated in the panels. The shaded area represents the la Poisson scatter 
around the mean. Observational data arc taken from Koposov et al. (2008) 
under the assumption of two different radial distributions of satellites, 
NFW-likc (solid circles with error bars) and isothermal (open circles). The 
arrows on error bars indicate that there is only one galaxy in that particular 
bin, and so the Poisson error is formally 100%. 



theless, this discrepancy is not as strong as that seen for 
the luminosity function and it does not allow us to reject 
the m w — 1 keV model. We should also point out that our 
determination of Mo. 3 can be partially affected by neglect- 
ing primordial thermal velocities that would create a core in 
the density profile (Strigari et al. 2006). The expected size 
of the core is 100 and 8 pc for the m w = 1, m w ~ 5 keV 
models respectively (Strigari et al. 2006), thus for m w — 5 
keV our approach for computing Mo.3could be not fully self- 
consisted. 



4 DISCUSSION AND CONCLUSIONS 

In this paper we compare the most recent observational data 
on the luminosity function of MW satellites with a suite of 
iV-body simulations combined with SAM techniques, with 
the aim of constraining a possible "warm" nature for DM 
particle. In particular, we want to make use of the new tight 
observational constraints on the faint-end, which have al- 
ready been proved to be of fundamental importance to test 
SAMs in the light of the problem of MW missing satellites 
(see M10) 

We perform a series of iV-body simulations with dif- 
ferent choices for the mass of the warm dark matter parti- 
cle (m m ), and we combine them with Semi- Analytical Mod- 
els (SAMs) for galaxy formation, using the best fit solution 
found in M10 for the parameters describing the SAMs. We 
then compare the resulting statistics of MW satellites with 
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Figure 5. Mass within 300 pc versus luminosity. Red dots show results 
from LCDM numerical simulations. The (blue) squares show results for 
WDM satellites for a m u mass of 5 kcV and 1 keV in the upper and lower 
panel respectively. Black points with error bars arc the observational results 
from Strigari ot al. (2008). 

both the predictions for a LCDM Universe and observational 
data. 

Our results show that we are indeed able to put a lower 
limit on m w (1 keV), which is in agreement with previ- 
ous determinations based on the galaxy power spectrum 
and Lyman-a forest observations (Viel et al. 2008, Boyarsky 
et al. 2009a). It is worth noting that our limits are less tight 
with respect to the previous determinations, as a number 
of caveats have to be taken into account. First, only the 
m w — 1 keV realization really fails on reproducing the MW 
satellites luminosity function (with our "standard" choice 
of SAM parameters). Using a m w value as small as 2 keV, 
we still predict a luminosity function which is almost indis- 
tinguishable from the concordance LCDM cosmology. This 
is due to the fact that the for m w > 1 keV, the small-scale 
power modifications do not affect considerably the substruc- 
tures responsible for the formation of the bulk of visible 
MW satellites, showing an intrinsic limit of this approach 
with the current observational data. Moreover, the faint- 
end of the LF, which provides us the strongest statistical 
constraints, is sensitive to the details of the description of 
baryonic physics included in the SAMs, and in particular to 
the strength of stellar feedback in low-mass halo (see M10). 
Therefore only models that drastically failed in reproducing 
the data (as m w — 1 keV) can be realistically ruled out. In 
the present work we only considered a very simple WDM 
model; it is worth notice that there are more complex and 
more physically motivated models discussed in the litera- 
ture (e.g. warm+cold dark matter, Boyarsky et al. 2009b 
or composite dark matter Khlopov 2006, Khlopov & Kou- 
varis 2008). These scenarios deserve special studies and can 
provide new, interesting hints on the nature of dark matter. 



This paper is mainly a ground test on the feasibility of 
our approach. New data coming from the most recent mul- 
tiwavelength surveys like the Panoramic Survey Telescope 
and Rapid Response System (Pan-STARRS jf| will provide 
a better determination of the MW satellite properties and 
then an optimal data-set to increase our knowledge about 
the nature of dark matter. 
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